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RANDOM MATRIX THEORY AND THE SIXTH PAINLEVE 

EQUATION 



Abstract. A feature of certain ensembles of random matrices is that the 
corresponding measure is invariant under conjugation by unitary matrices. 
Study of such ensembles realised by matrices with Gaussian entries leads to 
statistical quantities related to the eigenspectrum, such as the distribution of 
the largest eigenvalue, which can be expressed as multidimensional integrals 
or equivalently as determinants. These distributions are well known to be t- 
functions for Painleve systems, allowing for the former to be characterised as 
the solution of certain nonlinear equations. We consider the random matrix 
ensembles for which the nonlinear equation is the a form of Pvi . Known 
results are reviewed, as is their implication by way of series expansions for 
the distributions. New results are given for the boundary conditions in the 
neighbourhood of the fixed singularities at t = 0, 1, oo of crPvi displayed by a 
generalisation of the generating function for the distributions. The structure 
of these expansions is related to Jimbo's general expansions for the r-function 
of frPyi in the neighbourhood of its fixed singularities, and this theory is itself 
put in its context of the linear isomonodromy problem relating to Pvi . 



1.1. The (7- form of Painleve VI. Given a large sequence of an eigenvalue spec- 
trum, it is a simple matter to rcscalc so that the mean spacing between consecutive 
eigenvalues is unity, then to empirically determine the distribution function for the 
spacing. When these eigenspectra are the highly excited states of heavy nuclei, it 
is a celebrated result that the distribution function is well approximated by the 
functional form 



known as the Wigner surmise. It is in the exact computation of eigenvalues distri- 
butions for certain classical random matrix ensembles that Painleve transcendents 
make their appearance. These transcendents may relate to any of Pn to Pyi , de- 
pending on the random matrix ensemble under consideration, or its scaled limits. 
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1. Introduction 



(1.1) 




2 



P.J. FORRESTER AND N.S. WITTE 



As part of the centenary of the discovery of Painleve VI by Fuchs [11] , in this paper 
we will restrict attention to the random matrix ensembles relating to Pvi • 

Painleve VI conventionally refers to the four parameter second order nonlinear 
differential equation 

2 \y y-1 y-tj \t t-1 y-tj 

, y{y-l)iy^t) f /3t ^(t - l) Stit - I) 



as obtained by Fuchs [11]. However in random matrix theory this equation is never 
encountered directly. Rather what is encountered is the so-called Jimbo-Miwa- 
Okamoto cr-form of Painleve VI, 



(1.3) a'y, {t{t - lXi)'+(aVi [2avi - (2t - l)a'y,] + v^V2V^Vif ^ \{{(T'yi+vl) 

fc=i 



This is written in a form which displays a 1)4 root system symmetry in the param- 
eters 111, . . . , W4. When expanded out there is a common factor of cryj, which when 
cancelled out shows (1.3) to be a second order second degree nonlinear differential 
equation. The equations (1.2) and (1.3) are related by the Hamiltonian formulation 
of Pyi , due originally to Malmquist in 1922 [19]. 

In the Hamiltonian approach to the Painleve equations in general, one presents a 
Hamiltonian H(q,p,t; {vk}) where {vk} are parameters, such that after eliminating 
p in the Hamilton equations 

where the dash denotes derivatives with respect to t, the equation in q is the 
appropriate Painleve equation. The Hamiltonians can be systematically derived 
from the isomonodromy deformation theory associated with the Painleve equations 
[16], [22] (aspects of the isomonodromy deformation theory associated with Pyi is 
covered in Section 3 below). From such considerations, the Hamiltonian relating to 
Pyi was given by Okamoto [23] as 

(1.5) t{t-l)H^i = q{q-l){q-t)p' 

- [(f3 + Vi){q - l){q -t) + (W3 - Vi)q{q - t) - [vi -f V2)q{q - l)]p 

+ (W3 - Wi)(w3 - V2){q - t) 

where the the parameters wi, . . . , 1)4 are related to a, /?, 7, 5 in (1.2) according to 

(1-6)^ 1 1 1 

a = -(V1-V2Y, (i = --{v:i+Vif, ^=-{vi-V4)^, (5 = -(1 - (1 - -^2)^), 

and q satisfies (1.2). Note that iJyi is quadratic in p and thus according to the 
first of the Hamilton equations (1.4) p can be written as a rational function of q 
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and q', and thus in fact Hyi is a rational function in q and q'. According to the 
foUowing result, this particular rational function, augmented by the addition of a 
linear function in t, satisfies (1.3) [23]. 

Proposition 1.1. Define the auxiliary Hamiltonian 

(1.7) hYi{t) = t{t - l)i?VI + e2[-Wl, -V2,V3]t - ie2[-Wl, ~V2,V3,V4], 

where 

(1.8) ep[ai, . . . ,a,.] := ^ aj^aj^ ■ ■ ■ aj^. 

l<jl<---<jp<S 

This auxiliary Hamiltonian satisfies the a-form of Painleve VI (1.3). 

1.2. Historical Overview. There are certain ensembles oi N x N random ma- 
trices with complex Gaussian entries and invariance under conjugation by unitary 
matrices, which have their joint eigenvalue probability distribution function (p.d.f.) 
of the form 

1 ^ 

(1-9) ^n^2(x/) n i^k-x,)\ 

1=1 l<j<k<N 

where the weight function W2{x) is of one of the classical forms 

, Gaussian 

(1.10) ^2(2;) = x°-e~^ {x > 0), Laguerre ■ 

y'{l-x)'' (0<.T<1), Jacobi 

For example, let X be a n x N{n > N) rectangular matrix of complex Gaussians 
N[0, 1/a/(2)] + iN[0, l/\/(2)]. Then the matrix X^X has eigenvalue p.d.f. (1.9) 
with Laguerre weight x"~^e^^ {x > 0). The sixth Painleve equation relates to 

(1.9) with Jacobi weight; in particular to the probability that there are exactly n 
eigenvalues in the interval {t, 1) of that ensemble. This probability is in turn equal 
to the coefficient of (1 — ^)" in the expansion of 

(1.11) Ei{t-aMO-=^(^j^ )'^'''"' (i ^^l)'^''^ 

N 

1=1 l<j<k<N 

where C denotes the normalisation. In the case N = 1 (1.11) is an integral form 
of a particular 2F1 hypcrgcomctric function, which was related to the Painleve VI 
equation by Okamoto [23]. 

Let {Pj{x)}j=o,i^,,, denote the set of monic polynomials of degree j orthogonal 
with respect to the Jacobi weight x°'{l — x)''{0 < x < 1), which are given in terms 
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of the Jacobi polynomials pj"'''-' by 

(1.12) p,(.x) = (-l)^,!i>±i±l±l) p(-'')(l - 2x). 

Let {pj,Pj)2 Jo{Pj{^))'^^°'{^ ~ x)^dx, and define 

.o. i^J , ^ {w2{x)w2{y)Y'^ Pn{x)pN-i{v) - PN-l{x)pN{y) 

1.13) K^[x,y)^- ^ . 

It is well known, and easy to derive (see e.g. [4]), that with K"^ ,^ denoting the 



(1.15) t{t-l)-\ogEi 



integral operator on (t, 1) with kernel (1.13), 
(1.14) Eiit- a, b; = det(l ~ ^K'lj^^t^^^) 

It was in this form that 

d 

'di 

was first related to the solution of a nonlinear equation by Tracy and Widom [24]. 
The equation found was of third order. Subsequently Haine and Semengue [13] 
studied (1.11) itself in the case ^ = 1 and found a different third order equation 
for (1.15). Upon subtracting the two equations they obtained a second order sec- 
ond degree nonlinear equation which they identified as an example of the a-form of 
Painleve VI (1.3). The study of Tracy and Widom proceeded via functional proper- 
ties of quantities associated with the Frcdholm determinant (1.14), while Haine and 
Semengue used the theory of the KP hierarchy and Virasoro constraints satisfied 
by certain matrix integrals as introduced by Adler and van Moerbeke [1] . A third 
approach to the problem was initiated by Borodin and Deift [3]. They combined 
Riemann-Hilbert theory with the method of isomonodromic deformation of certain 
linear differential equations to obtain a characterisation of (1.15) which allows for 
immediate identification with the parameters in (1.3) (this is not the case with 
[13]). Explicitly it was shown that 

(1.16) a{t) = -t{t - l)j^\ogEi{\ - t- a, 5; ^ + viV2t + ^(--^1^2 + vi,Vi), 
satisfies (1.3) with 

a + b a + b a — b 

(1.17) Vi=V2=N^ —, W3 = , "^'4 = — 

Furthermore, it is required that as i — > 

(1.18) ^ logi?](,(l - t; a, b; - -CK'Jij{l -t,l-t) 

(1.19) =-CCAr(a,6)(l-<)', 
where 

n 90^ r ( h^- r(a + 6 + iv + i)r(6 + jv + i) 

(1.20) CN{a, b) _ r(iV)r(a + 7V)r(fe+l)r(5 + 2)' 
thus providing the boundary condition to be satisfied by (1.16). 
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According to (1.16) 



(1.21) E'^il^t;a,b;0^exp / ^ 



a{s) - V1V2S - -{-ViV2 + V^Vi) 



s) V 2 

From the characterisation of a{t) as a crPyi transcendent with a specific boundary 
condition, the power series sohition of (1.21) about t = 1 can be readily computed 
[4], 

(1.22) Eif{l~t;a,b;0 = l 

-C^(l-t)^- {l - ib + mN^ + 2ia^m-2-2b + ab)^^ _ ^ ^^^^ _ ^ 

, ^2 g^(a, b){N - 1){N + b + l)(iV + g - l)(jV + a + 6 + 1) _,^26+4 ^nn -fU 
+^ (& + 2)2(&2+46 + 3)2 i^+^l^ ^^t- 

Furthermore, one can anticipate from (1.14) that the leading term in 1 — i accom- 
panying the power ^'^ will be proportional to (1 — t)^^^^~ , as is consistent with 

(1.22) . 

It is well known (see e.g. [4]) that after the change of variables Xj = cos^ < 
9j < TT and with a, 5 = ±1/2 the eigenvalue p.d.f for the Jacobi ensemble as speci- 
fied by (1.9) and (1.10) becomes identical to the eigenvalue p.d.f. for matrices from 
the classical groups 0^{N), Sp{2N) chosen with Haar (uniform) measure. As a 
consequence (in an obvious notation) 

£;0-(2iv+i)((o^^).^) ^ i5^(cos2 0/2;C)| 

b=-l/2 

(1.23) i?o^(2^+i)((0,(/<);0 = E■^{cos'<|>/2■,0\.^~^/^ ■ 

fa = l/2 

Analogous to the expansion (1.22). the tiFvi evaluation (1.21) can then be used to 
deduce the expansions [4] 

(1.24) i?°-(^~+i)((0, x);i) = 1 - cx + ^^^cV 

36 

487^^ - 407V2 + 7 , - 5iV2 + 1 

cx^ H c^x^ 

3600 2025 

192iV6 - 3367V^ -I- 196iV2 - 31 , - 112N^ + 77 N"^ - li ^ ^ , 

H cx^ c^x^ + 0{x^), 

211680 198450 ^ 



(1.25) £;O-(^^+i)((0,x);C) = l-^^^cx^ 

(47V2 - 1)(127V2 _ 7) ^ g (4Ar2 _ i)(48Ar4 _ 72Ar2 + 31) ^ ^ 

3600 211680 ^ 

where c = 2N^/tt. 
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The method of [24] was adopted in [25] to relate 

1 / roc / 
(1.26) E'i'(,;n:i)-^=-{ I -( I )<ii, - (y -ij )dx 

1=1 ^ l ' l<j<k<N 



^ 1 



to crPvi • Exphcitly it was shown that 

(1.27) ais) = (1 + s')-^^\ogE'^y{s; a + N; 0, 
satisfies the equation 

(1.28) (1 + s^iay + 4(1 + s^){a'f - 8sa(a')' + 4.a\a' - a^) 

+ Sa^scra' + 4[iV(iV + 2a) - a^s^]{a'f = 0. 

As noted in [7], the relationship between (1.28) and (1.3) can be seen by changing 
variables 

(1.29) t^'-^^, a^^l^t)^ ^h{s), 
in the latter so that it reads 

4 

(1.30) h' ((1 + s^)h"f + 4 {h'{h - sh') - iviV2V3Vif + 4 ]^(/i' + vl) = 0. 

fe=i 

With 

(1.31) h = (T — a^s, Vl = —a, V2 = 0, = N + a, V4 = a, 

(1.30) reduces to (1.28). 

The interest in (1.26) in random matrix theory comes about by making a stere- 
ographic projection from the real line to the unit circle by the change of variables 
X = tan 6*72. With z — e*^ this shows 

^ 1 

1=1 ^ <■ ' l<j<k<N 

N 

= 2-^(^+2") n ii + ^if' n i^fc - • • • 

1=1 l<j<k<N 

In the case 77 = the measure on the right hand side of (1.32) corresponds to the 
eigenvalue p.d.f. for N x N random unitary matrices chosen with Haar (uniform) 
measure. For general 77 G N it corresponds to this same ensemble conditional so 
that there is an eigenvalue of degeneracy rj at 6 = tt. 

With £'^'^■'((01,02);^) denoting the generating function for the probability 
that the interval (01,02) contains exactly n eigenvalues, it follows that 

(1.33) E^^^\{0,2x);O ^cxp(^- h{cot (t))d(j^ , 
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where h{s) satisfies (1.30) with vi ~ V2 = fs = 0,V4 = N. Since for x — > 0, 
e'^^^\{0,2x);£,) ^ 1— ^xA^/tt, we seek the solution of (1.30) subject to the bound- 
ary condition h{s) ^ c,c :~ ^N/n. This allows the power series expansion 

(N'-mN^-2)i,N'-,) _ (iV-4)(iV-l)' 

52920 291600 ^ 

to be computed [4] . This expansion was first computed in [24] using the character- 
isation of i?^^^'' in terms of a third order nonlinear differential equation. 

We have given reference to three distinct approaches which relate (1.11) to non- 
linear differential equations with the Painleve property. There is a fourth approach, 
which is due to the present authors [5], [6] [7], and involves applying Okamoto's 
theory of the Hamiltonian systems approach to Pyi [23]. This approach has the 
advantage of allowing generalisations of the generating functions (1.11), (1-26) and 
e'^^'^\{4'1t 4'2)', to be related to crPvi • Consider for definitencss the latter. In 
[7] the more general quantity 

M., ..,.;..,...o - ^ /;j (/:-r /;j 

N 

1=1 i<j<k<N 

where t = e*"^, <j> £ [0, 2tt), ^* G C and the parameters uji, UJ2, n E C, oj ^ uji + iuj2, 
are restricted with 5R(2wi), SR(2^) > —1, N G Z>o. The independent variable 
t, whilst originally defined on the unit circle \t\ = 1 with a real angle 4>, can 
be considered as a complex variable which is analytically continued into the cut 
complex i-plane. The case 0-12 = /i = 0, cji = 77 of (1.35) gives the generating 
function for the probability of k eigenvalues in (0, 4>) for the ensemble specified by 
the right hand side of (1.32). Define 



(1.36) 

MN{a,b):= dxi... da;^TTz["-''/'|l + z,|"+^ TT \z, ~ Zk\' 

J -1/2 J -1/2 -,^,".7. 



^1/2 J-1/2 l<j<k<N 

(1.37) ^TT^r(a + 5 + , + l)r(, + 2) ^^^^2...,^,^ 

^ ^ Ti r(a + j + i)r(6 + j + i)' 

set ^* = 1 - (1 - ^)e""^ and denote by h{s) a solution of (1.30) with 



.,N, 



(1.38) 



Wl = -/i-CJi, V2 — il^2, V3 ~ N + IJ. + UJi, Vi = -fJ, + UJi. 
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It was shown in [7] that 



(1.39) AN{t■,u>l,u!2,^J,;C) ^ 



M^ifi + uj,fi + uj) 
Mn (0,0) 



X exp 



h{cot -) + LU2{N + LUi- /.() + (^i + ^i)^ cot - 



cW 



We remark that according to the theory in the sentence including (1.30) 

1 



(1.40) 



/i(cot-) 



-2icrvi( 



1 + 



A gcnerahsation of the Jacobi ensemble (l-H) may also be made in a similar 
fashion, by the introduction of the additional factor Yl^^i \s — xi\^ in the integrand. 
In the case ^ = it is shown in [7] that this generalised ensemble can be related to 
(1.35), and consequently that 

(1.41) avi{t) = (e2[-ui, -W2,W3] +Nfi)t- ^e2[-ui, -W2, 1's, ^4] - Nfi 



i(t-l)-logAjv(i;wi,W2,M;r =0), 



is a solution of (1.3) with the parameters 
(1.42) 



N- 



N - 



"2 



N 



V3 



N - 



V4 



2 2 ' 2 ' " 2 

Moreover a derivation of the transformation implied by (1.40) and (1.41) was given 
which tells us that the condition ^* = in (1.41) can be relaxed. 

The above mentioned generalisation of the Jacobi ensemble (1.11) has also re- 
vealed different Painleve connections to those found in [7]. In the mid 1990's, with 
^ = 0. this was studied from the viewpoint of orthogonal polynomial theory by Mag- 
nus [18]. In that study an auxiliary quantity occurring in the theory was shown to 
satisfy the Pyi equation (1.2) for appropriate parameters. We remark that further 
development of orthogonal polynomial theory in relation to (1.35) [9, 8] has been 
shown to relate to Backlund transformations in the Hamiltonian theory of Pyi , 
and to the so called discrete dPy equations. 

In the case 0^2 = A* = 0, t^^i = ^7 (1.35) is the generating hmction for the proba- 
bility that there are exactly n eigenvalues in (tt — 0, tt) for the eigenvalue p.d.f. on 
the right hand side of (1.35). In this case expanding (1.14) allow us to specify the 
s — > cx) boundary condition which must be satisfied by h{s). However, for generic 
parameters the boundary conditions were not given in [7]. Consideration of this 
latter problem is our concern for the remainder of the paper. A similar situation 
was recently rectified in [10] in relation to a class of multidimensional integral so- 
lutions of crPy [6] . The approach taken was to write the multidimensional integral 
as a determinant (involving the confiuent hypergeometric function), and to expand 
the entries of the determinant. An analogous strategy suffices in relation to (1.35), 
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now with the entries of the determinant given in terms of the Gauss hypergeometric 
function. This is done in Section 2. 

The form of the general expansion of a solution of the cPvi equation (1.3), or 
equivalently its r function about its fixed singularities t = 0, 1, c» has been given by 
Jimbo [15]. The exponents occurring therein are given in terms of the monodromy 
data associated with the isomonodromy deformation formulation of Pvi • Aspects 
of this theory and the results of Jimbo are revised in Section 3. 

One of the features of the general expansions is that they consist of two branches, 
whereas the expansion of (1.35) only exhibits a single branch. We show in Section 
4 that this is entirely consistent with the monodromy data associated with cryi; it 
is such that the coefficient in front of one of the branches vanishes identically. 



Consider the multidimensional integral (1.35). According to (1.40) (with the 
condition ^* = relaxed), we know how to relate its logarithmic derivative to a 
solution of the crPyi equation (1.3). However this property cannot be used to char- 
acterise the multidimensional integral unless an appropriate boundary condition is 
specified. As remarked in the second paragraph below (1.42) above, in [7] the inter- 
pretation of (1.35) for some special parameters as the generating function of a gap 
probability provided the boundary condition in those cases. However for general 
parameters no boundary condition was presented. 

Taking a different approach, namely the expansion of the elements in the deter- 
minant form of (1.35), the sought boundary conditions can be deduced as presented 
in the following result. 

Proposition 2.1. For generic values of fi,uj,LU the spectral average An has the 
following expansions. About t — subject to fi ~ Cu ^ Z we have 



2. Boundary Conditions 



(2.1) 



l + C 



2i simr (fi ~ uj) 



N N-1 




r(i -I- fc + ^ + ^)r(i + k- fi + uj) 



fc!r(2wi + fc + 1) 




X 



2i sin 7r(/i - u) + ^* e-^■'i^'~^) 

r(i + ^ + w)r(i + /i - (i')r(i + 2/i)r(7v - ^i + cd) 

T{N)T{N + n + i:})T{N + 2LJi)r2(2 - N + n ~ u) 
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About t = 1 subject to 2/i + 2uii ^ Z we have 

(2.2) A.^U , fe'r(2, + 2a., + fc + l) J N.j^-.) 

^ ' 't^i i-i- T(l + k + n + Lu)r(l + k + u + Q)] 2m + 2lji ' 



k=0 



^ (_l)JV+i / g-7rj(p-i.) ^ sin7r2^sin7r(/i + tj) 



sin7r(2/i + 2Li;i) \ 2i sin7r(2/i + 2^1) 

r(i + 2^l)T{l + 2^i)r(i + + ^)r(i + + _ 1+2^+2^, ' 

r2(2/i + 2wi + 2)r(2^ + 2cji + l)r(7V)r(-iV - 2/.t - 2wi) ^ ' 
And about t = oo subject to [i — uj <^'L we have 

(2.3) t-^^'A^^ Umn{f, + u;)+C ^ 

t— 00 1^ sin7r(/x — cj) \ 2i J 



JV-l 



fc!r(2wi + fc + 1) 



-XI + A: + /i + cj)r(l + fc - ^ + cj) 

^ r^2Ai7v(/.+<i)i 



N - /.i + uj t 

2i siii7r2^ + C*e~"^^ 
2i sin7r(^ + cu) + ^*e-'^*(A'+'-) 



^ r(i + 2^^)^{l + ^ + ^)r(i + ^ - ^)r(jv - /x + a;) i+Ar_^+^ | 
r(Ar)r(iv + M + w)r(Ar + 2LJi)r2(2 - iv + /i - w) j 

Any one of these boundary conditions suffices to uniquely define a solution to the 
ordinary differential equation (1.3) under the above generic restrictions on the pa- 
rameters. 

Proof. Using Heine's identity we can deduce from (1.35) that 
(2.4) A]s{t]UJl,uJ2,^i]C) 



= t-^f" det 



0<jM<N-l 



This is a determinant of a TocpUtz matrix whose symbol, or more specifically weight 
w{z), is one defined on the unit circle \z\ = 1 by 



(2.5) w{z) = r''z-^-"'(l + zf^^ (1 + tzf^' 



1, 6 e (-TT, - (f) 

ee(7r-0,7r) 



with Fourier components Wk such that w{z) = T^^=-oo'^i^^^ ■ Considering the ex- 
pansion about t = first we note that the Toeplitz matrix element can be evaluated 
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as 

r(2tJi + i) 



(2.6) t^u;„ = <il + e 



2i sin7r(7i + /i — cj) J r(l + n + /.j + uj)T{l — n — fj, + uj) 
X 2-Fi(— 2/i, —n — ^ — uj; 1 — n — iJ, + uj;t) 

e-'^'("+^-") r(2/i+l) 
2i sin7r(7T, + /.j — cj) r(l + n + /.j — t^)r(l — n + /i + iD) 

X t"+^-^2i^i(-2cji, ?i - - 1 + n + /i - 1:^; i), 

under the condition that n + /i — cj ^ Z, which is a form suitable for the development 
of an expansion about t = 0. In the course of the derivation we had to invoke tighter 
constraints on the parameters, namely \t\ = l,t 1 and 5R(2/^) > 0, K(2ti;i) > 0, 
but these can be relaxed by analytic continuation arguments. The structure is 
i^'uin = a„(t) + t"+''^'^fo„(t) where a„(<),6„(<) are analytic about t = 0. One can 
expand the Toeplitz determinant about i = retaining only the leading order terms 
from both the analytic and non-analytic contributions and the resulting formula is 

(2.7) t^^ANit) ^ det(a,_fc(0)+ta', fc(0)),,fc=o,...,JV-i 

+ (-l)^-i6_(jv-i) det(a,_fe+i(0)),, fc=o.....JV-2<''-^-^'^-'). 
Using the determinant identity [21] 

^9 8^ ^ r(d + fc-j) ^ ■, r(i + d-c) r(d-n + i + i) 

^'■^^ d^*l,r(c + fc-,)j„<^.,<„_,-ll^T(i + d-c-,) r(c + ,) ' 

or its slightly more general form 

'r(zfe + 6-jr 



(2.9) det 



r(2;fe - i) 



0<j,k<n-l 

n-1 



^ n (^>-^,)n(-i)'(-'') /'" n:)""" - 

0<j<fc<n-l j=0 ^ 

for an arbitrary sequence {^j}"^o^ one can evaluate the determinants appearing in 
(2.7). The resuh is (2.1). 

To develop the expansion about t = 1 the appropriate expression for the Toeplitz 
matrix element is 

(2.10) r(2, + 2., + l) 



T{l + n + fi + w)r(l -n + fi + Lu) 
X 2-F'i(~2a;i, n — fj. ~ Q; — 2/i — 2lui; I — t) 
7rt{n+,,-u) sin7r2^sin7r(?i + ^ + cj)\ r(l + 2/x)r(l + 2wi) 



2i sin7r(2^ + 2wi) J r(2 + 2/i + 2wi) 

(1 - f)i+2A'+2'^i2Fi (1 + 2^, n + 1 + + cj; 2 + 2^j + 2cJi; 1 - <), 
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subject to 1 + 2yU + 2u!i ^ Z. The structure is now t'^^'^Wn = a,i(l — + (1 — 
<)^+^^+^'^ifo„(l — i) where a„(l — 6„(1 — are analytic about t = I. Again one 
can expand the Toephtz determinant about t = 1 retaining only the leading order 
terms from both the analytic and non-analytic contributions and use the above 
determinant identities to arrive at (2.2). 

Lastly the expansion about t = oo \s computed using the Toeplitz matrix ele- 
ment in the form 

(2.11) t^Wn^-. -, sin7r(n + ^ + w) + r 



sin7r(n — /i + cij) \ 2i 

^ r(2^i + 1) 

r(l + n - /i + u})T{l -n + fi + oj) 

^-iri{n+tJ,+uj) / g-7ri2/i\ rf2/i + 1) 

' sin7r2^ + r ' ' v a- ; 



sin7r(n — p,-\- lo) \ 2i J T{\ + 7i + /i + — n + /i — w) 

X r+^+'^2i^i(-2wi, -n- n-uj;l-n + n^uj; \/t), 

valid for n - ^ -I- w ^ Z. The structure is i"^w„ = a„(l/t) 4- &„(l/t) where 

«n(l/i), b„{l/t) are analytic about t = oo. Repeating the procedure adopted about 
the other singularities one arrives at (2.3). □ 

Remark 2.1. We do not attempt to treat the degenerate cases where either fx — uj€ 
Z, 2fi + 2wi g Z or /i — ti> G Z here as this results in confluent logarithms and other 
technical difficulties. 



3. ISOMONODROMY DEFORMATION FORMULATION FOR Pyi 

In this section we describe the isomonodromic deformation system which char- 
acterises the general solution to the sixth Painleve equation and outline a solution 
to the direct monodromy problem, that is explicit formulae for the monodromy data 
associated with a particular solution of the a-form for the sixth Painleve equation. 
The isomonodromic deformation system associated with the sixth Painleve equa- 
tion is not uniquely determined and we shall see this arbitrariness arise in Section 
4 when applying the general theory to our random matrix theory. 

Following the conventions and notations of [15], [17] we consider the 2x2 linear 
matrix ODE for 5'(A;f;) with four regular singularities in the A-plane chosen to be 
= 0, t, 1, oo 



(3.1) 
(3.2) 
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It is taken that the residue matrices Av{t) satisfy 

Q 

(3.3) Ao + At+Ai = -A^^ = -^(73, a3:= 



1 

-1 



e c\z, 



and given the freedom to choose cither tvA^ or one of the eigenvalues we foUow the 
convention of Jimbo [15] 

1 



(3.4) 



trA^ = 0, det 



9t, J/ = 0,t,l,oo, 



defining the formal exponents of monodromy 9^. The residue matrices satisfy the 
Schlesinger system of equations 

[Ao,At] 



(3.5) 
(3.6) 
(3.7) 
(3.8) 



d_ 
di 
d_ 

di 
d 
di 



At 
Ai 



[Ao.At] , [AuA,] 



[At, A,] 
1-t ' 



0, 



as a consequence of the compatibility of (3.1) and (3.2). 
The T-function for Pyi is defined by 



(3.9) 

and the a -function 



dt 



logr = Tr 



Ao _A^ 

t t-1 



(3.10) 



C(t)=t(t-l)-log. + - 



which satisfies the second-order second degree differential equation 



dt 



(3.11) -C t(t-l)— C 



dt^ 



t-C-C 



dt \ dt 



dt 



1^' 



-rX + 7(0t + e^r -X + -i9t- 



dt 



dt 



(cf. (1.3)). 

Furthermore we suppose {Assumption 1) that the matrices A^, are diagonalis- 
able, i.e. that there exists nonsingular E SL{2,C) such that 



(3.12) 



R-^A,R, = -e^a^, 6, e C\Z. 
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In the neighbourhood of a regular singularity ^'(A) can be expanded locally as 



(3.13) * = E *"-(^ - i^)™+^"^a, 

m=0 



for ;/ = 0, i, 1 and for A = oo in the form 



(3.14) /+E*™ooA-™ A- 



For such local series to exist we have to assume {Assumption 2) that the eigenvalues 
of are distinct modulo the non-zero integers, i.e. that ±6i, ^ N. 

The matrices C^.v = 0,<, 1 arc the connection matrices and we are taking the 
local solutions (3.14) as our fundamental system of solutions, i.e. Coo = I- The 
monodromy matrices My[v ~ 0, t, 1, oo) are defined as 



(3.15) *lA=.+^e-. = *Ia=.+5 u = QX 1, oo, 

and are given in terms of the monodromy exponents and connection matrices by 

(3.16) = C-^e'''^"''''C^. 

The monodromy matrices together satisfy the cyclic relation 

(3.17) M^MiMtMo = /, 



according to the convention taken for the basis of loops displayed in Figure 1, which 
generate the fundamental group 7r(CP\{0, t, 1, oo}, Aq). There is arbitrariness in the 
monodromy data in the sense that the replacement C,y ^ D~^Ci, doesn't change the 
monodromy matrices provided that Di, commutes with the right-hand side of (3.12). 
This implies that is diagonal ii di, ^ 0. This arbitrariness will manifest itself 
in the appearance of an arbitrary complex number in the explicit parameterisation 
discussed latter. 
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Figure 1. Monodromy representation of the fundamental group 
for CP\{0,t, 1,00} 

The isomonodromic principle states that the monodromy data MD := {9,^, C^, M,j\v = 
0, 1, c«} are preserved under the deformations of t. The invariants of the mon- 
odromy data are defined 



(3.18) = 2cos7r6'^ := TrM^ /i e {0, 1, cx)} 

(3.19) 2cos7rCT^^ := TrM^il/^ e {0, t, 1}, 

in the sense that these do not contain any arbitrary constants. 

In [15] Jimbo states the following conditions under which his results apply 

(3.20) B^,euBi,e^iz, 

(3.21) < 5R((7ot) < 1, 

(3.22) e^±et±aM. eoo±ei±(JMi2z. 



When (Tot 7^ a parameterisation of the monodromy matrices was deduced by 
Jimbo and is given in Lemma 3.1. 



Lemma 3.1 (Jimbo[15]). Subject to the conditions (3.20), (3.21) and (3.22) the monodromy matrices can he parameterised in the 
following way 



(3.23) 







(3.25) CMtC-' = 

(3.26) CMqC 
where 



^' ' ^ ~ isin7r6loo \^2r"ie"''~ sin f(6loo - fi*! +cr) sin 1(6*00 - fi*! -fj) - costtct + e'^''^- costt^i j' g 

H 

_i _ 1 / e'^*'" cos 776*4 -COS 7r6lo -2se'^*'" sin f(6io + - cr) sin f (6*0 - 6lt + cr)\ > 



i sin TTcr \ 25-16-""" sin 1(6*0 + 6lt +cr) sin §(6*0 - - 0-) -e"''"" cos 776*4 + cos 776*0 /' « 



e''*'" cos 776*0 - cos 776*4 2s sin ^ (6*0 + 6*4 - cr) sin f (6*0 - 6*4 + cr) 



_ _L I O \j»JO /I l.'Q /I t.'i Dill 2\ (J^^C "^7 i^i^^ 2 

f(0o-^t-'T)sinf 



(3.27) C = 



-1 sin f (0OO - 01 + ct) sin f (0OO + Oi-a) J' 
Here r is an arbitrary non-zero complex number, and the short hand notation s := Sotjf = ""ot is used. 



•2. 



isin77(T \-2s-i sin f (6*0 - 6*4 - cr) sin ^(6*0 + 6*4 + cr) -e-""" cos 776*0 + cos 776*4 /' H 
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Proof. We begin by noting that (3.23) follows from our choice for the fundamen- 
tal system of solutions. To establish the other formula we require two prelimi- 
nary results. Let us make the abbreviations c^, :~ costtO^, :~ sin7r0^,Co. :~ 
cosTTCTot, S(T := sin TT (Tot, := -|- is^ and 



6(i9) sin-i^, 



(3.28) 

and note some trigonometric identities which will be useful in this and the ensuing 



proofs 




















(3.29) 


26(0oo - t 


'i + 




'oo + 


01 


+ cr) 


= ci - 


Ccc-Cff -\ 




(3.30) 


2&{e^ + i 


'i - 




'oo 


01 


— (t) 


= ci - 






(3.31) 






C^)S(^ 


'oo 


01 


— (t) 


= Ca - 


^^l^oo ~ 


Si Sqo , 


(3.32) 


2&{6^ + i 




(7)6{6 


'oo + 


01 


— a] 


= Ccr - 


ClCoc> -\ 


^l^OO J 


(3.33) 






C^)6(^ 


'oo + 


01 


— a] 


= -Coc 


+ ClC^ 


+ SiSa- 


(3.34) 


2&{6^ + i 




C^)6(^ 


'oo ~ 


01 


— a] 


= -Coc 


+ ClCcr 


- SlSa 


(3.35) 


26(00 - 


0t 


+ cr)6 


:oo + 


0t 


+ ^) 


= Ct~ 


CqCct + 


SoSa, 


(3.36) 


2e(6io + 


0t 


-a)& 


%- 


0t 


— a] 


= ct- 


CoCa - 




(3.37) 


26(00 - 


Ot 


+ cr)6 


:oo + 


0t 


— a] 


= -Co 


+ CtCa - 




(3.38) 


2S(6io + 


Ot 


+ cr)6 


%- 


0t 


— a] 


= -Co 


+ CtCa 





We also note that a general member M of SL(2,C), characterised by detA/ = 1 
and tiM = 2c, can be written as 



(3.39) 



M = 





where + s"^ — 1. and with arbitrary complex numbers b, b^^ ^ 0, x^^ ^ 0. One 
can parameterise c = cos7r0, s = sin7r0. The observation made in [2] that simplifies 
the derivation is that the Jimbo parameterisation is constructed so that 



(3.40) CMtMoC-^ = A := 

for a non-singular matrix C G GL(2,C) 



(3.41) CM^^M^^C-'^ 



Using (3.23) and the parameterisation (3.39) for Mi we can recast the above relation 
0jS cL homogeneous matrix equation for C, that is 




(3.42) 



C = C 



eoo(ci+a;) -eoo{si + ix)b\ 



si — IX 
be^c 



Ci — X 
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This has a solution if and only if ixsoo = Co- — ciCoo, which fixes x. Using the 
identities (3.29,3.30,3.31,3.32) the equation for the components of C can be written 

(3.43) -ieoo66(0oo + Oi + (j)Cn + 6(^00 - Oi - <j)Ci2 = 0, 

(3.44) eoo&6(0oo + O1- a)C2i + i&iBoo - ^i + o)C22 = 0, 

assuming Boo ± 0i ± cr ^ Z. If one defines an arbitrary complex constant r := ihtaa 
then this yields the solution (3.27) for C . We note that det C = — SooSo-- Employing 
the solution for x and the identities (3.31,3.32) we arrive at the formula (3.24) for 
M\. For the remaining monodromy matrices we utilise the representation (3.39) for 
CMqC~^ and CMtC~^ in the relation (3.40). From the (1, 1) and (2, 2) components 
of the resulting matrix equation we find 

(3.45) EooCo - ct = eoo2;o +a:t, 

(3.46) -(^cq + ct^ <^'^xq + xt, 

and one deduces that iSuXq — cqc^ — ct and is^xt ~ ctCa — cq. From the other 
components we find that x^ ~ x"^ ~ Cj — Cq and the formula for the ratio of the 
other undetermined constants is given by 

(3.47) T"- tir^7v 

bo &[9o-0t+cr) 

where we have utilised (3.35,3.36,3.37,3.38). If we define the remaining undeter- 
mined constant bt '■= —i^aSot then we recover the parameterisations (3.26), (3.25) 
forMo,A/t. □ 

A key identity is the following connection relation which relates sot,crot to an 
and a 01. 

Lemma 3.2 (Jimbo[15]). One of the connection relations is 

(3.48) 441 sin |(0o + OtT <^ot) sin |(0o ~dt± <Jot) 

7T TT 

X sin -(6*00 + CTot) sin -{9^ - Oi ± aot) 

_ g±7rjo-ot gjj-j^ TT(7ot COS TTfTtl — COS TrOf COS Tt9oo ~ COS TtBq COS tt6i ) 

± i sin TTdot COS ttctoi + cos irOt cos tt9i + cos tt^oo cos tt9q 

Proof. The proof of this has been detailed in Boalch [2] with a typographical correc- 
tion to the original formula in [15]. So we content ourselves with a brief summary of 
the steps involved. After noting (3.40) it was found that the monodromy invariants 
are more manageable in the forms 

(3.49) pti = tiiCM-^C-^A-\CMtC-^)), 

(3.50) poi - triCM^'C-'A-\CMoC-')), 
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because they were linear in 1, sot, Sot^- Upon taking the combination (3.49)+ecr(3.50) 
in order to eliminate the Sqj^ terms it was found that the resulting expression could 
be factorised through the use of the identities 

(3.51) &{e^+6i-a)&{6oo-0i-a) 

- e(0oo - 6*1 + Cr)6(6'oo + 6*1 + cr) = SocSa, 

(3.52) e^&{e^ + ei-a)e{e^-ei-a) 

- 6006(6*00 - dl+ a)&{9oo +61+ (t) = iSoo(ecrCoo - Cl), 

(3.53) e^&{e^+6i-G)&{6^-6i-a) 

- e^&iOoo -O1 + a)6{eoo +ei+<j)= isocici - e'^Coc), 

enabling a factor of 2sooSa to be cancelled out. This then yielded one of the desired 
formulae (3.48), whilst the other could be found from the other combination of 
(3.49,3.50). □ 

A consequence of this is a constraint on the monodromy invariants {pot,PtiiPoi} 
which is an algebraic variety defining a sub-manifold, the monodromy manifold, of 

Lemma 3.3 (Jimbo[15]). The monodromy manifold is given by 

(3.54) mipot,Pti,Poi) ■■= 

POtPtlPol +Plt+Ptl +-P0I - (PoPt + PlPoo)P0t - {PtPl +P0Poo)Ptl - {PoPl +PtPoo)P01 

+ PQ+Pt +Pi+ pIo + PaPtPiPoo -4 = 0. 

Proof. We multiply the upper and lower sign forms of the left-hand side of (3.48) 
to eliminate sot and then employ the identities (3.29,3.30,3.35,3.36) to replace the 
product of sines. Equating this to the corresponding product of the right-hand 
sides yields (3.54) as the only nontrivial factor. □ 

Remark 3.1. The above connection relation involves only the free parameters sot 
and the monodromy invariants (Tot, (Jti, coi. As we shall see immediately below the 
arbitrary parameters sot and (Tot appear in the expansion for the r-function about 
t = 0, and there exist analogous pairs about t = \,oo. Correspondingly there exist 
two other forms of the connection relation (3.48) involving either the parameters 
■Sti, soi and can be deduced directly from (3.48) by a simple substitution rule given 
at the end of Theorem 3.2. However both connection relations yield the same 
formula (3.54) for the variety defining the monodromy manifold. 

Now we come to the fundamental result for the expansion of the r-function in the 
neighbourhood of the fixed singularities of the sixth Painlcvc system at i = 0, 1, 00. 
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Theorem 3.1 (Jimbo[15]). Under the conditions (3.20,3.21,3.22) we have the ex- 
pansion of the T-function as t in the domain {t G C|0 < \t\ < s, |arg(i)| < (/)} 
for all e > and any (/) > 

(3.55) T{t) ^ ct^-'-ol-s^)/'^ 



8a2 



16(t2(1 + ct)2 



where a =/= and s are related to s through 

r2(i - (T)r(i + i(0o + dt + ct))t{i + i(-0o + 0t + ct)) 



(3.56) s^s 



r2(l + a)T{l + i(0o + et~ ct))T{1 + i(-0o + 0t^ ct)) 

r(i + ^(^oo + 01 + a))r(l + ^ (-gp, + 01 + a)) 
r(l + i(0oo + 01 - a))T{l + i(-0oo + 01 - a)) ' 



X 



and we employ the short-hand notation s = sot, s = .SQt and a — aot. The mon- 
odromy data defining the unique solution to the sixth Painleve system is {(JQt,sot}- 
Here C is an arbitrary constant. 

Proof. The details of this proof of this are given in Jimbo [15] and so we do not 
repeat them here. Also Guzzetti has laid out some of the intermediate steps in 
the appendix of his work on the elliptic representations of the general Painleve six 
equation [12]. □ 



The regular singularities x = 0,i, l,oo play equivalent roles and can be ex- 
changed under linear fractional or Mobius transformations. Consequently one can 
solve the connection problem very neatly and under the additional conditions 

(3.57) < 3fJ((7ti),5R(aoi) < 1, 

(3.58) 01 ± 0t ± an , 0oo ± 0o ± <Jti i 2Z, 

(3.59) 00 ± 6li ± aoi , 0oo ± 6»t ± aoi ^ 2Z, 



derive expansions about t = 1, cx). 
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Theorem 3.2 (Jimbo[15]). Under the conditions (3.20,3.57,3.58) we have the ex- 
pansion of the T-function as t ^ 1 

(3.60) r(i) ~ C(l - t)(-?i-«?-''?)/4 

- 16.?,(l + -n)^ 

and as t oo 

(3.61) r(t) ^ 

f //)2 _ /i2 _ ^2 \ //)2 _ /)2 _ ^2 



X {1 + 



i - g"oi)(^o ~ ^1 ~ "'Ol) ,-1 



9^ - (g«-aoi)^][g^-(gi-aoi)^],_i_.„, 



16a2,(l + aoi)2 



- (^t +C^Ol)^][^0 - (^1 +gOl)^] .-l+gQi , n/U|-2(l-3?(<T0l))N 



16fTo^i(l-aoi)2 

Here Sti,soi are found by making the following substitutions in (3.56), (3.48) re- 
spectively 

(3.62) s ^ sti, s^sti, 9o ^ 01, a ^ ati, aa ^ aot, 

(3.63) s^Sqi, s^sqi, 9q^6oo, cr ^ o-qi, o-qi^CToi, 

with 
(3.64) 

COSTTfJoi = — COS TTfTot ~ 2 COSTTCToi COS TTfTtl + 2(cOS TT^Q COS TT^t + COS TT^qo COSTT^l). 

T/ie monodromy data defining the unique solution to the sixth Painleve system is 
either {ati,sti} or {<toi,Soi}. 

4. Monodromy Data for the Spectrum Singularity Ensemble 

The precise relationship between the spectrum singularity average AN{t;) and 
the isomonodromy theory of the sixth Painleve system is given by the following re- 
sult. Its validity relies on the conjecture that the expansions of Jimbo given in The- 
orems 3.1 and 3.2 remain valid upon relaxation of the constraints (3.20,3.21,3.22) 
and (3.57,3.58,3.59), provided the former are well defined (i.e. do not then diverge). 



Proposition 4.1. For the spectrum singularity ensemble the associated isomon- 
odromic system is not unique but the monodromy data for any of these systems 
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falls into three generic cases. An example of each case is given below in cases (A), 
(B) and (C). The formal monodromy exponents can be taken to belong to either of 
three sets 

(4.1) Case(^) : 6*0 = -/i - t^, 6*4 = ^^ + 2^1, 6*1 = TV + 2^i, eoo = ~t^-^, 

(4.2) Casc(S) : Oq^h-Oj, 9t ^ N, Oi = N + 2fi + 2uji, 0oo = tJ.-i^, 

(4.3) Case(C) : Oq = -2uji, 9t= N + fi + u, 6*1 = iV + /i + cj, 600 = 2^1. 

The monodromy invariants for either case are 

(4.4) CTot = iV — /.J + cj, tJti = 2/i + 2cji, (Toi = iV - /i + cj. 

In the case (A) the monodromy coefficients are 

2isin7r(^ — Q) 



(4.5) sot = 1 



(4.6) sti = l + r 

(4.7) SOI = - 



^-iiiip-u,) sin7r(2^ + 2wi) 



2i sin7r2yLtsin7r(/Lt + w) ' 



- 1 + e-^'^v- 

All monodromy matrices are lower triangular 



(4.8) A/o = 



(4.9) M, 

(4.10) Ml 



^ mo e'^'(''+'^)^ 

^g7ri(JV+2i^i) Q 



* I mt e-'^^f^+^-i) i ' 



'g7ri(Ar+2M) Q 

-7ri(Ar+2/i) 



mi e 



where 



(4.11) mo = — 



2i r sin7r2cji sin7r(/i + w) siii7r2/isin7r(/i + w) 



sin7r(/^ — tj) [ Sot r 

(4 12) mt = 2ii(-l)^sin7r2u;i f sm7r(/i + cj) _^ sin7r2^ 



sin 7r{fi — Lo) [ sot 
(4.13) ^, = -Mzi)!^£!!iI^e-(-+-) 



For Case (B) the monodromy coejficients are 

2i sin7r(/i — w) 



(4.14) sot = 1 



sin7r2a;i sin7r(u + ij-') sin7r2/isin7r(u + u;) e "'(^ 
(4.15) Sti : — \ ^ = ^-^7- ^^^^ — r-^+C 



sin7r(2^ + 2a;i) sin7r(2/i + 2a;i) 2i 

(4-16) SOI = ^,_^^^4.»p • 

One of monodromy matrices is proportional to the identity, the others are full 
(4.17) 

^ i / e^'^*'^"") cos7r(/^ — Lij) — cos7r(2^ + 2wi) 2r sin7r(/^ + Lij) sin7r2/x 

° sin7r(/i " w) I sin7r(^ + w) sin7r2wi -e"'^^^'^) cos 7r(/x - w) + cos 7r(2^ + 2ti;i) 

(4.18) Mt = i-lfl, 

(4.19) 

^_ i(-l)^ /e-"(^-") cos7r(2^ + 2011) -cos 7r(/x-w) 2re-''''('"-'^) sin7r(^ + (I>) sin7r2^ \ 

^ sin7r(^-a;)i |e''*(^"'^) sin7r(^ + tj) siii7r2tJi C0S7r(^ - cj) - e''*^''"'^) cos7r(2^ + 2cji) i 
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- 1 + e'i'^^M 

All monodromy matrices are upper triangular 



4.20) 
4.21) 
4.22) 



For Case ( C) the monodromy coefficients are 

2isin7r(/i — Q) 



sti- 



' ^*p-T^l(^i-Q) ' 

sin7r2wi sin7r(/L< + (D) sin7r2/isin7r(/.t + tj) 



sin7r(2/x + 2uji) 



soi = 



sin7r(2/i + 2wi) 



~Yi ' 



4.23) 
4.24) 

4.25) 

where 

4.26) mo 

4.27) mt 



Mt 
Ah 



y e 





mo \ 

— 7ri2cJi 1 ' 

mt \ 
mi \ 



2i 



sin7r(/i — uj) 
2i(-l)^sin7r(^ + w) 



{— sin7r2/isin7r(/i + cj)sot + sin7r2u;i sin7r(/i + uj)r} 



|sin7r2^e""(''"'^'sot - sin7r(/i + (I>)r| 



sin7r(/i — ti>) 
4.28) mi = 2i sin 7r(^ + (I.)e""(^+2^V. 

Proof. Comparison of the two differential equations for the ^-function, (1.3) and 
3.11), imply that in general 



4.29) 



{vi,V2,V3,V4} = - < 



£3(^0 + ^1) 
[64(^0-^1) 



with Cj = ±l,j = 1,2,3,4 and eie2e3e4 = 1- Using either set of parameters, (1.42) 
from the JUE correspondence or (1.38) from the CylJE correspondence, we find 
that the monodromy exponents can be given by one of three sets 

{N + 2n, N + 2lui, fi + oj, fi + UJ 
N,N + 2fi + 2uji,fi~Lj,^i-Lj 
N + fi + uj,N + fi + Q, -2^, 2wi 

modulo permutations of the monodromy exponents and an even number of sign 
reversals. This is a manifestation of the non-uniqueness of the isomonodromic 
system for our problem. For definiteness we choose one example of the three cases. 
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namely the cases (A), (B) and (C) given in (4.1,4.2,4.3). We note some simple 
identities which do not depend on the choice of the permutation or the sign 

(4.31) vl + vl + vl + vj = ^{el + el + 61 + el), 

(4.32) v,v,V3V, = j^iel - el){9l - 91). 

and in particular the following products which apply equally to cases (A), (B) and 
(C) 

(4.33) 

i^l - ODidl ^ 9^) ^ {N - fi + uj){N + ti-uj){N + 3iu + uj){N + iu + 2oj + uj), 
(4.34) 

{9l - 9l){el -9l) = {N-fi + u;){N + ^ - i:^){N + 3/i + oj){N + iJ + 2oj + oj), 

(4.35) {e'l - 9l){el - 9l) = {2N + 2fi + 2uj{){2^i - 2uj{){2^i + 2uj{){Cj - uj). 

If we make a comparison of the r-functions themselves for the JUE correspondence, 
(3.10) and (1.41) we find, at the level of the r-functions 

(4.36) Ajv(t;) = 

Now applying the t Q expansions for T{t;9), namely (3.55), with those of 
AN{t;), (2.1), we first note that the exponent of the t-prefactors must be consistent 
and this implies 



(4.37) = 1(91 + 9^ + 91 + 91) + 2e,[v^n = (^ u™^)^ = {N - ^, + u)^ 

i=i 

This result applies to all the cases and for definiteness we make the choice of sign 
outlined in (4.4) (the other choice of sign is essentially equivalent). Turning to the 
leading analytic term of order t in both expansions we find that its coefficient is 
given precisely by 

2n{fj, + Lj)N 



(4.38) 



for all three cases upon employing our solution for ctqj. Next we make a comparison 
of the non-analytic terms in (2.1) and (3.55) and it is here that we have to treat 
the cases separately. However it is generally true that in the classical situation of 
the finite rank random matrix ensemble only one of the non-analytic terms is ever 
present, the other being switched off through the following mechanism. Taking 
case (A) first our task is to show how the coefficient of the non-analytic term t^+'^»t 
vanishes and to match the remaining coefficient with that of the application. The 
identification of dot immediately implies dot = 9q + 9t. Therefore 

(4.39) [ti+'^'«] oc „ }. = 0. 
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From a comparison of the coefficients of the remaining non-analytic term we find 
precise agreement and this enables us to determine the solution for the monodromy 
coefficient sot in (4.5). For case (B) we note that aot =^ 0t ^ &o and thus 

(4.40) [ti+'"'»] oc — ^ 0. 

Similarly we find agreement for the coefficient of the surviving non-analytic term 
and that enables us to fix the monodromy coefficient as in (4.14). In case (C) we 
have (Tot ~ 0i — 0oo and thus 

(4.41) oc 3-473 = 0. 

Again we find agreement for the coefficient of the t^+'^o* term and conclude that the 
monodromy coefficient is given by (4.20). For each of the three cases we observe that 
the Jimbo parameterisation of the monodromy fails (see condition (3.22) however a 
meaningful result emerges so we conjecture that the Theorem still holds with these 
relaxations. 

Now we make a comparison of the expansions at t = 1, namely (2.2) and (3.60). 
Examination of the algebraic prefactor using the relation (4.36) leads us to conclude 
cr^i = (2/i -I- 20^1 )^ in all three cases and we choose the positive sign, as in (4.4). 
Employing this solution we compute the coefficient of the analytic term is, in all 
three cases, 

(4 42) 

which is are entirely consistent with that in (2.2). To examine the non-analytic 
terms we take the three cases separately again. For case (A) we see that the 
coefficient of the (1 — t)^""^" vanishes because an = —Oq — O^o and 



(4.43) [il-tf-"*^] oc n \ ^ =0. 



The coefficients of (1 — t)^+'^" now agree precisely provided we have the solution 
(4.6) for the monodromy coefficient sti- For case (B) we have the relation ati — 
9i — 6t and see that 

(4.44) [(1-t)^-""] oc = 0. 

Again the coefficients of (1 — f)^+'^" agree and the solution (4.15) for sn follows. 
In case (C) we see that crti = 0oo — Oq and this ensures 

(4.45) [(1 -i)!-'""] oc TT-L-z =0. 

In this case wc also find the coefficients of the remaining non-analytic terms are 
precisely consistent, leading us to deduce the solution (4.21) for sti- 
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It remains to make a comparison of the expansions at i = oo, namely (2.3) and 
(3.61). Using (4.36) we see the algebraic prefactor implies that ctoi ~ ^ l^ + ^Yi 
and we choose the positive sign for the exponent. Using this value for ctoi we 
compute that the coefficient for the term in all three cases is 

(4.46) 2,i,^^)N 

which is consistent with (2.3). To treat the non-analytic terms we take the cases 
separately. For case (A) we note that ctoi = 6t + 9oo and this implies 

(4.47) [r^-""^] oc „ }„ = 0. 

The coefficient of the i-i+'^oi term is found to be in agreement with that of (2.3) 
if we take the solution (4.7) for sqi- For case (B) the relation is (Tqi = St — Ooo and 
this in turn implies 

(4.48) [r^-""^] (X = 0. 

Again exact agreement is found for the other coefficient provided that (4.16) holds. 
Lastly in case (C) we have the same relation as above and the absence of the t~^~"'°^ 
term. Examination of the coefficients of t~i+°'oi then lead us to the solution (4.22). 

Now we come to consideration of the connection relation (3.48) for t = and its 
two equivalent forms for t = 1, oo with respect to our solutions for the monodromy 
data. We compute that the three connection relations of either sign decouple into 
a left-hand side and a right-hand side which vanish separately for all the cases (A) , 
(B) and (C). The left-hand sides for the t = connection relation vanish because 
da + 9t - <7ot = and 6*00 + 6*1 + (Jot = 2iV for case (A), 6*0-614-1- aot = and 
6*0 + 6*4 + (Jot = 2N for case (B), and 0oo - Oi + (iQt = and 6*0 + 6*4 + aot = 2iV 
for case (C). Similar reasoning applies to the connection relations at t = 1 and 
t = 00. The right-hand sides of the relations vanish identically for both signs with 
the evaluations of aot, Cti, aoi as given in (4.4). 

To conclude we compute the monodromy matrices for the three cases and note 
that case (A), case (B) and case (C) yield the classical monodromy structure of 
lower triangular matrices, full matrices with one being a signed multiple of the 
identity, and upper triangular matrices respectively. □ 

Remark 4.1. Our results are consistent with the findings of Mazzocco [20] which 
state that the classical non-algebraic solutions for Pyi have either reducible mon- 
odromy groups (cases (A) and (C)) or at least one monodromy matrix is equal to 
±/, that is the monodromy group is 1-smaller (case (B)). Both these cases cover 
the situation of a one parameter (N) family of classical solutions. 



Remark 4.2. We observe that the exponents aot, <^ti, foi are not free boundary con- 
ditions for classical solutions but are fixed by certain combinations of the formal 
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monodromy exponents. Related to this phenomenon is that aU the connection rela- 
tions decouple so that the coefRcients of the monodromy coefficients sot, sti, soi all 
vanish and thus cannot be determined from these relations. There is a geometrical 
picture of the classical solutions, which was discussed in relation to Painleve II by 
Its and Kapaev [14]. The classical solutions of Pyi define singular points in the 
monodromy manifold which are characterised by SUl = and 

PtlPoi + 2pot - PoPt - PlPoo = 0, 

PotPoi + 2pti - ptPi - PoPoo = 0, 

PotPti + 2poi - PaPi - PtPoo = 0. 
We verify that these relations are satisfied for the cases (A), (B) and (C). 
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